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ABSTRACT 
A theoretical analysis for the characteristic undamped natural 


frequencies and mode shapes of three dimensional piping systems with 


general topology is presented. The development makes use of matrix 
methodology including transfer matrices, dynamic stiffness matrices, and 
certain topological matrices derived from the theory of linear graphs. 


A distributed mass model is employed throughout. 
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CHAPTER I 
Introduction 
1.1 General Remarks 
The dynamic behaviour of piping systems is of general interest in the 
design of such systems. The objective of this thesis is to present an 
analytical method of determining the characteristic frequencies and mode 
shapes of mechanical vibration of three dimensional piping systems possess- 
ing general topology. 
Although no digital computer program is incorporated in this thesis, 
the matrix methods employed and the iteration solution suggested not only 
lend themselves to a computer solution, but they make a large capacity, 


high speed digital computer mandatory. 


1.2 Scope of Work Presented 

A method for determining, analytically, the undamped natural fre- 
quencies and mode shapes of three dimensional piping systems with general 
topology is presented herein. Previous investigators in this field have 
presented techniques which have been successfully employed for analysis of 
systems contained entirely in a plane and with "tree like" topology, that 
is, elastic branches have been permitted, but involved meshes or loops have 
been excluded. Also, previous mechanical models have not permitted the 
analysis of circular, elastic arcs with distributed mass. A method of 
obtaining the transfer matrix for a circular arc with distributed mass is 
presented in Appendix A. The details of the method of obtaining topological 
matrices which permit analysis of general topology is presented in Chapter 
2, Section 2.7. Dynamic stiffness matrices are defined in Chapter 2, Sec- 
tion 2.4. The details of the necessary geometric analysis are presented in 


Chapter 2, Section 2.5 and 2.6. 





1.3 Notation. 


[ ] matrix 

[sv ] state vector 

SV, i th element of a state vector 

[2] 4 transfer matrix from node a to node b 

[D], deflection partition of a state vector 

(f], force partition of a state vector used with transfer matrix 


sign convention 
[Z,,],ete six by six partition of a transfer matrix 
sales dynamic stiffness matrix between node a and node b 
[S,,],ete six by six partition of a dynamic stiffness matrix 
[F]. force partition of a state vector, associated with 


dynamic stiffness matrix sign convention 


(Py. ob deflection matrix of the terminal nodes of a primary path 
ae force matrix of the terminal nodes of a primary path 
m2 2 system coordinates of a node 
L vector length 
XK, Sad direction cosines 
eC unit vector 
o central included angle 
72) radius of curvature 
AO vector from point A to point O 
[G] node rotation matrix 
[TLG] state vector transformation matrix from a local 


reference system to the global reference system 





[TGL] 


[SGL] 


[SLG] 


state vector transformation matrix from the global 
reference system to a local reference system 

primary path deflection matrix transformation 

matrix from global reference system to a local reference 
system 

primary path force matrix transformation matrix from 
local coordinates to global coordinates 

network 

node incidence matrix 

primitive stiffness matrix 

node stiffness matrix 

frequency determinant 

mode shape matrix 

node force matrix 

node deflection matrix 

eigenvalue of an array of coefficients of 


differential equations 





CHAPTER II 
Definition of Properties and Operations 
2.1 Problem Nomenclature and General Procedure 

The solution for the characteristic frequencies of a piping system 
with general topology requires a form of network analysis. In the analysis 
which follows, a given piece of pipe will be referred to as an element. The 
ends of an element will .be incident on nodes. If only two elements are in- 
cident on a node, that node will be referred to as a trivial node. When 
more than two elements are incident on a node, that node will be referred to 
as a primary node. When an element terminates at a foundation or some other 
boundary of the system, that node, upon which only one element is incident, 
will be referred to as a boundary node. Those elements connecting two non- 
trival nodes together will form a primary path when only trivial nodes are 
in the path connecting the two non-trivial nodes. 

The necessary solid geometry analysis utilizes conventional vector 
notation and involves many different frames of reference. Right handed 
triads will be developed which will be associated with each end of each 
element and will generally be referred to as local triads or local refer- 
ence systems. Such a triad will be so oriented that the x coordinate will 
be tangent to the centroidal axis of the element involved. One triad will 
be designated a global triad or global reference system and the node at 
which the global reference system is located will be designated the global 
origin. 

The general procedure is as follows. Beginning at a primary or bound- 
ary node and proceeding along a primary path, certain geometrical properties 
will be developed for the first element encountered which will permic. comput- 
ing a transfer matrix. The next element will then be analyzed for its 
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geometrical properties, its transfer matrix will be computed, a node rota- 
tion matrix, which is a function of the two triads associated with the ends 
of the elements incident at the trivial node involved, will be computed, and 
the product of these three matrices will be formed so as to eliminate the 
trivial node and leave a single transfer matrix for the first two elements. 
The next element is then analyzed in a similar fashion, a new transfer matrix 
is formed, and so on until a single transfer matrix has been computed for 
the entire primary path. This transfer matrix will then be transformed into 
a dynamic stiffness matrix which will contain the characteristic dynamic 
properties of the primary path analyzed. The dynamic stiffness matrix will 
then be transformed to the global reference system. This procedure will be 
repeated for all primary paths. 

Once dynamic stiffness matrices have been calculated for each primary 
path, they will all be assembled into a single diagonal matrix known as the 
primitive stiffness matrix. A topological matrix known as a node incidence 
matrix will then be formed. By premultiplying by the transpose of the node 
incidence matrix and postmultiplying by the node incidence matrix itself, 
the primitive qaaeenees matrix may be transformed into a new matrix called 
the node stiffness matrix. This matrix will relate the forces at the 
boundary nodes to the deflections at the non-trivial nodes of the network. 
By applying a qualitative knowledge of the constraints at the boundary nodes, 
a determinant known as the frequency determinant may be obtained from the 
node stiffness matrix. The frequency determinant will vanish for the char- 
acteristic frequencies of ae network. Essentially, by a process of assum- 
ing values of frequency and determining the corresponding values of the 


frequency determinant, its zeros may be located and refined. 





2.2 State Vector 

The state vector at a eiven node, a, of an elastic system is the 
collection or group of quantities whose sense and magnitude completely 
describe the instantaneous displacement, both rectilinear and angular, of 
that point from its quiescent position, as well as the corresponding recti- 
linear and angular forces in the member at the same point at the same in- 
stant, that is, those forces which, if a cutting plane were passed through 
the point at the instant of interest, would have to be applied to each cut 
face to prevent relative motion between them. For the purposes of this 
discussion, the state vector will be considered to be a column matrix. 

Although the systems to be studied are three dimensional, the individual 
elements which make up the systems will be planar, that is, only straight 
elements or circular arcs are permitted. For such planar elements, it is 
possible to consider an in-plane and an out-of-plane state vector. Each of 
these state vectors would have six elements, three pertaining to force and 
three to deflections. Although this notion has some utility when discuss- 
ing systems contained entirely in a plane or when deriving the distributed 
mass transfer matrices (to be discussed later) of the individual planar 
elements, it is useless for describing a three dimensional system. There- 
fore, all future references to a state vector, symbolized by [SV], will be 
the complete three dimensional state vector containing twelve elements, six 
of which describe displacement, and the other six describing force. These 


elements are organized and defined as follows: 


r¢ 





SV 


[sv] = SV 


SV 


A local 


1, displacement in an x direction 


2, displacement in a y direction 
3, displacement in a z direction 
4,.rotation about an x direction 
5, rotation about a y direction 
6, rotation about a z direction 2-1 
7, force in an x direction 

8, force in a y direction 

9, force in a z direction 

10, moment about an x direction 
ll, moment about a y direction 
12, moment about a z direction 


x y z coordinate system and the elements of a state vector are 


arranged as in Fig. 2.1-l. 





Fig. 2.1-1. 





The x axis of a local coordinate system is tangent to the centroidal 
axis. The y axis is normal to the x axis and directed in a convenient manner. 
The z axis is normal to both the x and y axis a6 as to form a right handed 
triad. As an example, the y axis of a curved element may be fixed in the 
plane of the curve and directed towards the center of curvature. The z 
axis would then be normal to the plane of the curve. A more complete dis- 
cussion of coordinates is contained in a general discussion of solid geo- 


metry which follows. 


2.3 Transfer Matrix 

The derivation of transfer matrices is treated in great detail in 
Reference 2 and the method employed to develop the transfer matrix for a 
circular arc with distributed mass is discussed in detail in Appendix A. 
Accordingly, only the briefest discussion will be given here, largely for 
the purpose of refreshing the reader's prior knowledge of this subject. 
For the purposes of this discussion, a transfer matrix is defined to be a 
twelve by twelve matrix which relates a state vector at one node and with 
a given orientation to another state vector at another node in a common pri- 
mary path and generally having a different orientation from the first. Let 
[Z] represent such a transfer matrix. This relationship is described by 


the following matrix equation. 
[Sv], = [2] 5 ° [sv], 2-2 


In this equation, [Sv]. represents a column state vector at node a, 
[SV], represents another column state vector at node b, and [Z] represents 
the transfer matrix between node a and node b and is a frequency dependent 


property of the element or elements connecting node a and node b. 





A dot written between two matrices indicates nothing more than ordinary 
matrix multiplication which could, of course, be indicated by writing the 
matrices in close juxtaposition. Here the dot is used to avoid confusion 
which might result from typewritten typography. 

Consider the following system containing two elements designated ab 


and bc whose transfer matrices are known as [Z] 4 and [Z} 0. 


Fig. 2.-2-l. 
Given the state vector at node a, [Sv] > the state vector at node b, 
[Sv], is given by eq. 2-2. 


Furthermore, [SV] | is given as follows: 
[SV]. = [2], ; [Sv], 2=3 


It is apparent that the intermediate state vector at node b may be 


eleminated and the following equation written: 


[SV] = [2]... : [2] : [SV] | 2-4 
or 
[sv] | = [Zz]. : [SV]. 2-5 
where [Zz]. ~ [Z},. : [2], 2-6 


2.3a Transfer Matrix Sign Convention 
A sign convention, unique to transfer matrix operation, must be estab- 


lished. The force elements of [Sv], are considered to be applied to element 
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ab at node a by node a. They are positive when their sense is in the posi- 
tive x, y, or z direction of the triad associated with element ab and located 
at node a. The angular forces, or moments, are positive according to a right 
hand rule about the x, y, or z directions. The force elements of [Sv], as 
computed by eq. 2-2, are those forces applied by element ab to node b and 
have sense and orientation according to the triad associated with element 

ab and located at node b. The force elements of [Sv], as given in eq. 2-3, 
are considered to be applied to element bc at node b by node b. Their sense 
and orientation are specified according to the triad located at node b and 
associated with element bc. In a like manner, the force elements of [Sv]. 
are applied to node c by element bc and have sense and orientation accord- 
ing to the triad associated with element bc at node c. 

It is apparent that there will generally be two local reference systems 
at each trivial node and that they will generally not have the same orienta- 
tion. This problem in geometry is handled by a rotation matrix which has 
been omitted from eqs. 2-4 and 2-6 in the interests of clarity, but is in- 


cluded in a discussion of geometry which follows. (cf. Sect. 2.8) 


2.4 Dynamic Stiffness Matrix 

A given state vector may be partitioned into two smaller column matrices, 
one containing the first six elements pertaining to deflection and referred 
to as a displacement vector, [D], and a second column matrix containing the 
last six elements pertaining to forces and referred to as a force vector, 
[f£]. 

Consider two nodes, a and b, connected by a primary path. The path 
connecting nodes a and b need not be a unique one. If there exists more 


than one path, this discussion applies to any given one. Let the two state 
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vectors, [SV]. and [svi be partitioned into the associated displacement 
and force vectors, [D] > (f] > [D],» and [f],- Let the appropriate transfer 
matrix, [2] be partitioned into four six by six matrices as indicated be- 


low. 


| D . | Zu | a D 2-7 
fe b Zo) | Zar ae a | 
Equation 2-7 may be expressed in an equivalent form as the following two 


equations: 


[D],. = [Zy ] e [D], + [Zz] ° [tf]. 2-7a 


(f], = (22,] . (Dl, + [221 - (£1, 2-7b 


Let the following six by six matrices be defined as indicated. 


ee) = tz) . tare 2-8 
[Syp] =-12, 17" 
(Spy) = [254] - (yg) + Wy)» (244) 
[Soo] = [Zo] - [241° 
We may now write the following. 
-1€], = (8,,] - D1, + [S,,] . (I, 2-9 
if], = [S.,] : [D], + [S.5] : [D], 2-10 


2.4a Dynamic Stiffness Matrix Sign Convention 

Equation 2-9 computes the negative of the forces at node a applied by 
node a to element ab. If it is clearly understood that the negatives of 
these forces are the positive forces applied by member ab to node a, the 
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negative sign in equation 2-9 may be eliminated and both (f]. and (f], are 
then considered to be forces applied by element ab to nodes a and b oriented 
according to triads located at nodes a and b. It is emphasized that this 
constitutes a unique sign convention-to be associated with stiffness matrices 
versus transfer matrices. In the succeeding development, force column 
matrices will be indicated by capital letters instead of lower case letters 
when they are associated with the stiffness matrix sign convention. 

The two force vectors may now be combined into a single twelve element 
column matrix containing the forces applied by the element ab to the nodes 
at a and b. In like fashion, the two displacement vectors may be combined. 
If the four six by six matrices are also combined to form a twelve by twelve 


matrix, the following equation may be written. 
| | _ 


ine, M12 | Pla 2-11 
aly PPT a2) lb 

which we write as . 
Fl, ,7 (sl, Pl, 2-1la 


The twelve by twelve matrix, [S], will be referred to as the dynamic 
stiffness matrix. The dynamic stiffness matrix, like the a matrix 
from which it was derived, is a function of frequency and the relative 
displacement and orientation of the state vectors at a and b and also depends 
on the element or elements connecting a to b by a given primary path. One 
can think of eq. 2-lla as relating the distortion of the member, [D]. pb? to 


the forces, [F]. b? producing them. 
9 
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2.5 Local Solid Geometry 


Analysis of a three dimensional piping system requires, for each 


element, a knowledge of the length, radius of curvature (for circular arcs), 


relative position and relative orientation. Such information may be readily 


obtained if the positions of the ends of each element are specified in a 
convenient basic reference system of cartesion coordinates. Such coordine- 
ates will be referred to as system coordinates and might be measured from 
a corner of a room or any other convenient reference point. The three 
coordinates will be referred to as X, Y, and Z. It is convenient also to 
develop a number of local coordinate systems for use in calculating trans- 
fer matrices and for other purposes. It is the purpose of this Section to 
accomplish this development. 

Consider a single straight element whose X, Y, and Z coordinates are 
known. Let one end be designated an "'a'' end and the other a “'b'' end. The 
"a" end may be considered to be the end with which a known state vector is 

TPT 


associated. Let the system coordinates of the “a” end be designated Xo» 


Yo and Za In a like manner, the “b'' end coordinates will be Xi Yh? and 


Zee Three local direction numbers may be defined as follows. 
X ob - x - Xo 2-12 
ob = Xb - Y. 2-13 
Z ob = Zo - Zs 2-14 


These three numbers are the direction numbers of a vector directed 
form "a" to "b'' and of length, L, equal to that of the straight element. 


L is given by the following equation. 





Three direction cosines Ky, Px and Xx may be obtained in the follow- 


ing manner. 


x ab 2-16 
L 

fhe, 2-17 
L 

te tab 2-18 
L 


These three direction cosines are the direction cosines of a unit 
vector tangent to the centroidal axis of the straight element. The direc- 
tion of this unit vector, designated a is defined to be the x direction 
of a unit triad fixed in space at the "a'' end of the element. For a straight 
element, the y direction of the unit triad is defined so as to be parallel 
to the XY plane of the system coordinates as well as being normal to the 
local x direction. Given a unit vector, Ey: directed in such a manner, 


the following vector equations may be written: 


ey aig |! 2-19 
e .e #=9Q 9-20 
x be 
where 3 a i + R.3 + ok and a = are ce &,3 


For F vector with a zero Z component, equations 2-20 and 2-22 may be 


written in scalar form as follows: 


Oa ey 2-19 
y mS % 
mK % XK, + B*R = 0 2-20a 


These two equations define the magnitudes of the direction cosines at ee 


as follows: 


2 2-21 














Cy moO 2-23 


These equations will be valid except when e, has zero Y component. 


In such a case, the following equations will be used: 


KX =0 2-2la 
y 
= ] 2-22a 
fy 
XY =0 2-23a 
y 


The z direction of the local triad at "a'' may now be defined by the 

following vector equation: 
e. = eX e 2-24 

The local triad thus established for a straight element specifies the 
orientation of the two state vectors associated with the ends of the element. 

The second of the two types of elements permitted in a three dimension- 
al piping system is a circular arc. In general, a piping system would be 
fabricated so that the ends of a circular arc would be tangent to the ele- 
ments that they are connected with. Consider such an arc with ends design- 
ated "'a"' and ''b"' and with the system coordinates of the two ends designated 


X , Ye Z Xs Ye and Z, in the same manner as the preceding straight 


a a b 


element. The element attached to the arc at the "a" end will have a triad 
associated with it whose orientation has been previously determined. The 
unit vector, ec.» of this triad will be tangent to the arc at the "a" end. 
Let this vector now be designated e,. Let the vector C connect the end 


points "a" and “b". The direction numbers of C, C Cys and C, are given 


1’ 5 
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by the following equations. 


Cy - X - Xo 2-25 
C, = X - a 2-26 
C, = Za = Z. 2-2/7 


The length of C is given as follows: 


Z Z 2 
L = Cy a C. ts C, 2-28 


The direction cosines of C are as follows; 


Os ao I279 
E 

f= C, J=30 
L 

Y Gp 2-31 
- 


These three direction cosines define a second unit vector e. contained 
in the plane of the arc just as is eee A unit vector, e» normal to the 


plane of the curve can now be defined by the following vector equation. 


eo= “eX 2-32 
n => 

t ‘ c 
A fourth unit vector, ee again contained in the plane of the curve, 


but directed towards the center of curvature, can now be defined as follows: 
= 2x = 2-33 


Designating the direction cosiness of @ as > &..» and ve the 
central included angle of the arc, @ , is given by the following equation. 
(Refer to Fig. 2.5-1l). 


@= 2 arcsin (e .e ) 2-34 
r Cc 
The radius of curvature, A, is given as follows: 


PP = L/2 2-35 
oe ° ae 
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Fig. 225-1 


The three unit vectors, eu e” and eo form a local triad with origin 


at "a'' which specifies the orientation of the state vector associated with 
the "a" end. A similar triad fulfilling the same function at the "b" end 
may be obtained as follows. 


— 


AO = AB + BO 2-36 

pe, = Le. +B 4 2-37 

e z@ - Le 2-38 
rb GC ic 

2 =6@ xX @é 2-39 


where eb? e , and é@ forms the local triad at the "b" end. 


n rb 
A combination of two straight elements and a circular arc, such as 
Fig. 2.5-2, where one is given the system coordinates of three points, 
"al, "b", and "c'', and the radius of curvature of an arc to be fitted in 


the included angle is also permitted. Since the required analysis is a 


combination of the previous two, it will not be presented here. 








2.5a Node Rotation Matrix, [G]. 
It is now possible to discuss the geometry requirements associated with 
more than two triads located at one common node, but with different orienta- 


tions. 





Fig. 2.5-3a Fig. 2.5-3b 

Conemder two elements, ab and bc, connecting three nodes, a, b, andc, 
with node b a common node. If the centroidal axes of the two elements 
are not collinear atnode b, as indicated in Fig. 2.5-3a, the triad associated 
with element ab and located at node b will not have the same orientation as 
the triad associated with element be and also located at node b. In order 
to perform the multiplication indicated in eqs. 23, 2-4, and 2-6, it is 
necessary to perform a rotation transformation upon the state vector [sv], 
as computed by eq. 2-2, before utilizing it in eq. 2-3 to compute [Sv]. 
The circular are transfer matrix developed in Appendix A performs the re- 
quired rotation associated with a curved element as well as the translation 
associated with all transfer matrices. 

Consider some element of [sv], as computed by eq. 2-2, oriented in 
the x direction of the triad associated with element ab at node b. [It will 
have three components when expressed in the primed coordinates of the triad 
at node b associated with element be. Since the triads concerned are de- 
fined by unit vectors expressed in system coordinates, the component of 
any given x directed element in the x' direction will be the product of the 


magnitude of the x directed element and the dot product, e .ée. 


q 
x x 
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As an example, let the first three deflection elements of [SV], 
oriented with element ab have magnitudes svl, sv2, and sv3. Let the desired 
magnitudes of the corresponding elements of [sv], oriented with element bc 
be designated svl', sv2§ and sv3' The latter three values are computed by 


the following equations. 


svl'=é'. @ . svl +@'. @ . sv2+6'.@ . sv3 2-40 
x x x° Zz 
sv2' =e'.@é .svl+éeé'.é . sv2+6'.é@ . sv 2-41 
y x y y 2z 3 
sv3' = e6'.@ . svl+éeé'.e . sv2+e'.e . sv3 2-42 
zx Zz zZ Zz 
or 
sv1L' @'.é e'.é e'.é svl 
x" x x y x° z 
sv2' |= j|e'.e e'.e e'.e sv2 2-43 
y Xx Viaaeay. y 2 
sv3'! e'.é e'.e e'.e sv3 
zx zy zz 
Or 
svl1' svl 
sv2° |=! fex] . sv2 2-43a 
sv3' sv3 
where 
e'.e e'.e e'.e 
x’ x x y Zz 
[GX] = 2-44 
e'.e e'.e e'.e 
y XxX y y y 2 
e'.e e'.e e'.e 
Zz x Zz y z° 2] . 


Since the four sets of x, y, and z directed elements of a given state 
vector all require a similar transformation, it is apparent that a single 
12 x 12 matrix will perform the desired transformation for all elements. 
Such a matrix will have the three by three matrix, [GX], repeated four 


times on the major diagonal as follows: 
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[G] = 2-45 


o ©" wo ex) 
The indicated elements of [G] are three by three matrices themselves. 
The bracket notation for [GX] has been omitted for clarity. This practice 
will be continued in the succeeding development. 
Now the following complete transformation can be made 
[sv]'= [G]) . [sv] 2-46 
Equations 2-4 and 2-6 can now be written without omitting the nec- 


essary rotation matrix: 


[sv]_ = [2] 21, 


(Z],, = (21, - (Gh, + (2), 2-68 


° [SV] | 2-4a 


2.6 Global Solid Geometry 

In order to make use of certain features of linear graph theory which 
apply to mechanical network analysis, it is necessary to express all directed 
quantities in a single common reference system. Such a reference system will 
be called a global reference system. A discussion of network analysis of 
mechanical systems follows in Section 2.7. 

The location and orientation of a global reference system may be chosen 
artibrarily. Although it could be anywhere, for the purposes of this dis- 
cussion, the global reference system will have the same orientation and 
position as the first local triad to be constructed. 

The derivation of the necessary transformation matrices and a more 
complete discussion are presented in Ref. 3. They are simply defined and 
exhibited here. 

Let [F] be some force column matrix expressed in a local reference 
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system and [F]' be the same force matrix expressed in the global refer- 
ence system. [F] and [F]' will be related by a transformation matrix as 
follows: 
(F]' = [TOA] . [F] 2-47 
The matrix [TOA] may be partitioned into four three by three matrixes 


as indicated below. 








TOA,, | TOA, 
| TOA5, | TOA, 5 
TOA, 5 is a zero matrix. 
TOA), is obtained precisely like [GX] of Sect. 2.5 with the exception 


that the primed coordinate system is now the global reference system. 


TOA = TOA 2-48 


TOA,, is equal to TOA, premultiplied by the 


following translation matrix [L]. 


0 -L, L, 
L = L, 0 -L, 2-49 
Ty 1 4 


The elements of L are defined as follows. 
Let X, Y, and Z be system coordinates of the local reference system 
position. Let m2 “P and an be the system coordinates of the global 


reference system position. Then 


Lo & - xX 2-50 
1 g 

L. =Y-Y I= 5 
2 g 

Vetere o *Z 2-52 
3 g 


An inspection of [TOA] shows that the inverse is given as follows: 
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feoeitoay = | 1 ° 2-53 
TAO, , TAO, , 
Where 
T 
[TAO, |] [TAO,,] = [TOA, |] 2-54 
and 
[TAO,.] = [TAO,.] . [L]° 2-55 
21 1) aoa 
The inverse transformation relationship can now be written as 
| 
[F] = [TAO] . [F] 2-56 


A similar transformation exists for displacements. 
Let [D] be some displacement matrix and let [D]‘' be the same displacement 
matrix expressed in a global reference system. The following equation may 


be written: 


[D]}' = [rao]~ . [Dp] 2-57 
where a 
TOA TOA 
T 11 2l 
[TAO] = re rng 2-58 


The inverse relationship is also true. 


[Dp] = [TOA]. . [Dp]! 2-59 


TAO TAO.. | 
Ty: i a 
[TOA] = li TAO, 2-60 


If the forces and deflections associated with a transfer matrix equa- 


where 


tion are transformed to a global form, a similar transformation must be made 
on the transfer matrix. A complete state vector is transformed as follows; 


TAO! 0 


[sv], = P . [SV]. 2-61 
0 TOA, 
oa [sv]' = [TLG], . [sv], 2-61a 


Ze 





TOA’ 0 
[sv] = - [sv]° 
0) TAO, 2-62 
= ’ - 
Or [SV]. [TCL] , , [Sv] | 2-62a 


These transformations operate on the transfer matrix as follows: 


[sv], = [2], - [sv], 2-63 
[SV], = [TLG], . [SV],= [TLG], . [2] _,- (T6L], . [Sv] ° 2-64 
or 
[sv], = (21), - [sv] 2 2-64a 
where 
[2]! = (TL¢], . [2],, - [TeL], 2-65 


A similar transformation exists for the dynamic 
stiffness matrix. The deflections are transformed as follows: 


TOA 0 : 2-66 


T 
[Dp] , = - ; {Diy 
a,b 0 TOA a,b 
bs 
Or [D]. b = [SCL], 4 ; Was 2-66a 
TOA. 0 
[Fle tb = : Ns 2-67 
0 TOA, 
or By > = OEE Le, . eee 2-67a 
Flo - [SEG] | ° web : yt) Os ° (0) 3b 2-68 
or Eb = ae , [D] ) ob 2-68a 
where 
pan = ISLE]. , IS}, > ; [SCL], 2-69 


In the preceding development, one may think of the transformation 


symbols as follows: 
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[TLG], local to global transformation matrix for transfer matrices 
[TGL], global to local transformation matrix for transfer matrices 
[SLG], local to global transformation matrix for stiffness matrices 


([SGL], global to local transformation matrix for stiffness matrices 


2.7 Application of Graph Theory and Topological Matrices 

The present treatment of the problem of Srerion of a mechanical net- 
work with general topology utilizes certain features of the theory of Linear 
Graphs. This discussion will be limited to those features and will be 
developed in terms of mechanical parameters rather than abstract mathe- 
matical ones. 

Given all the necessary quantities expressed in a global reference 
system, the topological solution which follows makes use of three conditions. 
First, the forces, applied by the primary path to the nodes upon which they 
are incident, are related to the end deflections of the paths by dynamic 
stiffness matrices. Second the nodes of the network must be in equilibrium, 
that is, the forces applied by the incident elements to a node must sum to 
zero for a primary node, or sum to the force which the node applies to the 
boundary in the case of a boundary node. Third, continuity requirements must 
be met at the nodes, that is, the deflections of the ends of the primary paths 
incident on a node are identically equal to each other and are eaual to the 
deflection of the node. 

Consider some network designated G. Let the number of nodes in G be M. 
Let the number of elements connecting the nodes be N/2. The N force column 
matrices associated with the N ends of the N/2 elements may be designated as 
[FE,], [FE,], [FE,],---, [FE]. The deflection matrices associated with the 


ends of the N/2 elements may be designated as [DE,]; [DE,], [DE,],---, (DE, j. 
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As defined by the stiffness matrices sign convention, the forces are 

positive when applied to the nodes. Now consider some node, m, of G, con- 
nected to the rest of G by n elements. In order that the node m be in 
equilibrium, the n forces associated with the elements incident on node 

m must sum to zero or to the force applied by the boundary to node m. The 
N-n forces that are not applied to node m do not enter into the summation. 
Let the boundary forces, which may be zero, be designated [FN J. There are M 
summations to be made, one for each node. A given summation may be indicated 


as follows; 


— C — 
[FN] a Jae ; [FE 5] 2-70 


The N coefficients [a']_ j are either a six by six identity matrix 
b 


or a six by six zero matrix according to the following rule. 


[a = j = I if end j is incident on node m 2-71 
[a] j = 0 if end j is not incident on node m 2-/1a 
b 
When the M summations have been made, the coefficients [aé] . form 


m,J 


a matrix designated [a]® which has 6M rows and 6N columns. Now, if the N 
force matrices are assembled in a single column matrix, [FE], and a second 
column matrix [FN], is formed which contains all the boundary forces, both 
zero and non-zero, then [FN] is related by equilibirum requirements to [FE] 
as follows: 

[FN] = [A]° . [FE] 2:72 

Next we concern ourselves with geometrical continuity at the same 

node m with n elements incident on it. The n deflections of the n element 
ends at m must be identical. Also, the deflection of the node itself must 
be identical to each of the n element deflections. Let the node deflec- 


tion be [DN J. One of these identical relationships may be expressed as 
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the following summation. 


M 


(DE) = 4, (kK), + DN 2-73 


The M coefficients oe are either a six by six identity matrix or a 

six by six zero matrix according to the following rule. 

a = [I] if end j is incident on node m 2-74 
oe = [0] if end j is not incident on node m 2~-74a 

There are N summations to be made. The coefficients Tele form a 
matrix which is [K]. Inspection of the rules for forming [K] reveals that 
[K] is the transpose of a. Accordingly, it will be referred to as [A] 
henceforth. 

If the N deflection matrices associated with element ends are assembled 
into one column matrix, [DE], and the M deflection matrices associated with 
nodes are assembled into one 'column matrix, [DN], then [DE] and [DN] are 
related by continuity requirements as follows: 

[DE] = [A] . [DN] 2-75 

As developed in Section 2-4, eq. 2-12a, the forces and deflections 
associated with a given element, e, having ends a, b, are related as 
follows: 

Sida = [SoJaip Daas 2-76 
There are N/2 such equations which can be written for G. Let them 


be written as follows: 


~(D, ] 


Fila om was l’a,b 276.1 


Fo) a.b im [Soda (Doda 2-76.2 

[F4), ob [Saab “(Data b 2-76.3 

Fy j2ta,p * [Soyo]. y «(Dyyo]. .2-76.N/2 
: N/2° a,b “zee COC 
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The N/2 equations may be written in more compact form as follows: 
[FE] = [SP] . [DE] 2-77 
In eq. 2-77, the ordering of [FE] and [DE] follows a certain method. 
The forces occur in pairs for each element. This imposes a limitation on 
how the topological array [A] is made up. Although the order in which nodes 
occur is not restricted, it will be found most convenient to take the bound- 
ary nodes first and then the remainder in any order. The order in which 
the forces [FE] have been numbered is governed by the manner in which the 
matrix [SP] is assembled. 
If eq. 2-75 is substituted into eq. 2-77 for [DE], the following is 
obtained. 
(re) = [SP] . [A] [DN] 2-78 
If eq. 2-78 is premultiplied by aie and eq. 2-72 is used, the follow- 
ing is obtained. 
[FE] = [A]" . [SP] . [A] . [DN] 2-79 
which we will write as 
[FN] = [SN] . [DN] 2-79a 
where 
[SN] = [A])" . [SP] . [A] 2-80 
The matrix [SP] will be referred to as the primitive stiffness matrix. 


The matrix [SN] will be referred to as the node stiffness matrix. 


2.8 Boundary Conditions and the Frequency Determinant. 

A knowledge of the boundary constraints provides certain qualitative 
knowledge of [FN] and [DN] which in turn makes it possible to extract from 
[SN] an array, designated [SN], the determinant of which will vanish for 
certain values of frequency which are natural frequencies of the physical 


system. 
ZT 








Fig. 2.8-1 


If the problem consisted of the simple beam of Fig. 2.8-1, the two 
nodes a and b would both be boundary nodes. The deflections at b would be 
generally non-zero while the forces would be zero. At node a, the deflec- 
tions would be zero while the forces would be non-zero. Depending on the 
boundary conditions at a given node, there may be any combination of zero 
and non-zero forces and deflections. However, there will always be the same 
number of zero quantities, that is, for each zero force there will be a non- 
zero deflection and vice versa. Reference 1 gives a more complete dis- 
cussion of such boundary conditions. 

The relationship [FN] = [SN] . [DN], as developed in Section 2.7, is a 
global equation while the given qualitative knowledge of [FNj and [DN] is 
in local coordinates. Just as forces and deflections were expressed in a 
global reference system so as to facilitate the desired topological opera- 
tions, they must now be expressed in local coordinates so as to make use of 
the given qualitative knowledge of boundary conditions. However, only those 
six element matrices containing both zero and non-zero elements in local 
coordinates need be expressed in local coordinates. If [FN] and [DN] have 


not been assembled so that those six element matrices just mentioned are ad- 


jacent and in low numbered positions, let [FN], [DN], and [SN] be reordered as 


necessary so that the qualitative information is so positioned. If the 


number of nodes is M and the number of those having non-zero [FN matrices 
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is p, then [FN] may be partitioned into two column matrices, the first 
containing p six element matrices and the second containing (m-p) six element 
matrices, of which, all elements are zero in global or in local coordinates. 
[DN] may be partitioned in like manner. The given qualitative knowledge 
of the first p six element matrices is not available in global coordinates, 
but the second will contain m-p six element matrices, which will be non-zero 
in either global or local coordinates. 

Let [FN], [DN], and [SN] be partitioned in this manner. Then the node 


stiffness equation* may be written as follows: 


Bai] ig]. Pe 2-81 
0 mile | S159 Des 

[F = [sn,,] : iD! + [sn,,] . bin 2-81la 
[0] = (sn,,]- (1 + [sn,,] . Dp! 2-81b 


The forces and deflections associated with the p boundary nodes have 
the same local position and orientation as the forces and deflections associ- 
ated with the element end incident on a given boundary node. Therefore, the 
inverse of the P translation and rotation matrices, [TOA] and [TAO] ©, which 
transformed the boundary node element end forces and igeneeetens to a global 
reference system may be retrieved and utilized to transform [F and [D,! to 
local coordinates. Let this be done and let the appropriate transformation 


matrices be arranged in two large arrays as indicated below: 


TAO 0 0 


1 
0 TAO, 0 = [TAOP] 2-82 
0 0 TAO 
P 


kof, Eq. 2-79a. 
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TOA 0 0 


l 
t Ct 
0 TOA, 0 = [TOAP] 2-83 
0 0 Toa’ 
p 
Then we can write 
F = [TAOP|] . [F 2-84 
! p! local ! 2 
and 
t 
[D,] = [TOAP]” . 1D.) yoca1 2-85 


Premultiplying eq. 2-8la by [TAOP] gives the following result. 


[TAOP] . [FJ = iv loos [TAOP] . [sn,,] : e + 


[TAOP] . [sny,] - [D,_)] 2-86 


Substituting eq. 2-85 into eq. 2-81b and eq. 2-86 for [DJ gives the 


following two equations. 


ft E 
iF! loca [TAOP] . [sn,,] - [TOAP] . [D,} local + 
[TAOP] . [sn, 5] : [Dp! 2-87 
C 
[Oo] = [sn,,] - [TOAP]~ , ID, Ts + [sn,,] : [D-p! 2-88 


Let the following four arrays be defined as indicated. 


[sn,,]" = [TAoP] . [sn] . [TOAP]* 2-89 
[sn,.]' = [TAOP] . [sn,,] 2-90 
[sn,,]' = (sn,,2 . [ToaP] © 2-91 
[sn,o]' = [sn] 2-92 


Now equations 2-8la and 2-81b may be written as follows. 
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= (omy . = 2-93 


oT 


m-p 
where el and fe) are now in local coordinates. 

There are 6p elements in [F! and also in ic Just as [SN] was re- 
ordered so that boundary node forces and deflections were in low numbered 
positions in [FN] and [DN], [SN]' may be reordered so that the non-zero 
elements of [Fo] are in low numbered positions and the zero elements of 
[D1] are in corresponding positions. For m nodes, [SN]' would be a 6m by 
6m matrix. If the number of non-zero elements of LF were q, then [SN]' 
could be partitioned as follows: 


F] = [sn,,J]°' - [0] + [sno ]ue | [D__¢] 7-94 


[Oy Sesiisn,.1°* . [0] + [sn,,]"' - Dee 2-95 


21 m-q 


For the boundary conditions given above, [sn,,]"' is q by q, [sn.,]'' 


12 
is m-q by m-q. Inspection 


is q by m-q, [sn,,]'' is m-q by q, and [sn,,] 


7g | 


of [sn reveals that it must be singular if the non-zero elements of 


92] +8 


[DN] may have any given value. The determinant of [sn will be referred 


ool"! 
to as the frequency determinant, |SNI . By evaluating |SNI at various 
frequencies over a given range, those values of frequency which yield a 

zero value for |SN| that is the eigenvalues or characteristic frequencies, 
may be determined. 

The reordering procedure indicated in the preceding is most appropriate 
for computer procedures. If one were doing such a manipulation with pencil 
and paper, one would naturally strike out rows of [SN]' corresponding to 
non-zero elements of [FI and strike out columns of [SN]' corresponding to 


zero elements of (DJ. The resulting array would be [sn We may now 


ool": 


adopt a more convenient notation for [sn,,]'', that is [SN]. 
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2.9 Mode Shapes 
If we presume that [SN] is evaluated at a characteristic frequency and 
that [SN] is only “singly singular", that is, its rank is one less than 
its order*, we have 
O = [SN] . [DM] 2-96 
Let [DM] have g elements. One, say B,» can always be chosen such that 
the other g-l elements of [DM] can be expressed in terms of g.° [DM] can 


then be normalized so that 


[DM] . [DM] = 1 2-97 

The matrix, [DM], thus obtained is unique, that is, independent of our 
choice of g.° 

[DM] is the mode shape of the primary nodes. Trivial nodes of interest 
may be included in [DM] by treating them as primary nodes in all the pre- 
ceding analysis. 

An optional method of obtaining the mode shape of the trivial nodes is 
as follows: 

Reconstruct [DN] including the elements of [DM] where appropriate and 
carrying out the necessary local to global transformations. Obtain [DE] by 
pre-multiplying [DN] by [A] as follows: 

[DE] = [A] . [DN] 2-98 
The matrix [DE] may then be partitioned into the twelve element column 


matrices associated with the primary paths. 


“otherwise, we have the coalescence of two roots which would complicate 
subsequent discussion unnecessarily, the situation being essentially no 
different than the case of any vibrating system having identical character- 
istic frequencies, corresponding to which the modes are not unique. 
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For any given primary path, say between primary nodes a and b, we may 


then write 


[ms (S] ie 2-99 


Solving eq. 2-99 for [F]. and (FJ, provides sufficient knowlédge to 
assemble a complete state vector at either node a or node b. By choosing 
node b, we may avoid any confusion concerning sign conventions. We may 


then write [SV], as follows: 


(ee =D) 2-100 


AES 


In eq. 2-100, we use lower case f indicating forces and sign conven- 
tions appropriate for transfer matrix use. We may then evaluate the 
transfer matrices, at the appropriate eigenvalues, between the primary nodes 


and proceed to obtain the deflections associated with all nodes. 
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CHAPTER III 
Discussion 
3.1 Network Analogies 

As has been previously stated, the solution for characteristic frequen-~ 
cies of three dimensional piping systems with general topology requires a 
form of network analysis. Once all vector like quantities had been exe 
pressed in a single common global reference system, as suggested by Ref. 3, 
the writer found the methods of neceoeh analysis developed by Fenves and 
Branin, Seshu and Reed, and Kron to be applicable to the subject analysis. 

Kron, in Ref. 5, describes the equivalent network of a vibrating beam 
as a six wire, six phase transmission line. Such a model does not permit 
the direct application of the development of Ref. 3. However, a similar 
technique was found in Ref. 4 which was directly applicable to a mechanical 
network when the properties of such a network were described in terms of the 
dynamic stiffness matrix suggested by Pestel in Ref. 1.* 

It should be noted that the mechanical transfer matrix was found to be 
analogous to the electrical transfer matrix of h parameters. The dynamic 
stiffness matrix was found to be analogous to the electrical short circuit 
admittance matrix. The connection tensors of Ref. 4 are node incidence mat- 
rices. And finally, force and current, as well as voltage and deflection, 


are analogous. 


3.2 Practical Limitations 
In principle, there are no limitations to the size of a mechanical system 
that can be treated in this manner. 


*cf. page 150, para. 3 of Ref. 1 
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In reality, there are two. First, any given digital computer, though possess- 
ing total recall, has a limited memory capacity. The method developed in 

this thesis requires the assembly of very large arrays. A network composed 

of 20 non-trivial nodes and 30 primary paths requires storage capacity for 
57,600 quantities for just two of the matrices involved, no matter how clever 
the programmer.* Although the subsystem technique developed in App. B would 
relieve this requirement somewhat, very large high speed computers are still 
mandatory computational tools. The second limitation is that of the size of 
the numbers involved. A very large system will require the ability to 

compute very large or very small numbers. Thus, double, or even triple, 


precision may be required in order that precise answers may be obtained. 


* 

For the example above, Allee [SP] would be dimensioned 120 by 360 

and would contain 43,200 elements. [A]© . [SP] . [A] would be dimensioned 
120 by 120 and would contain 14,400 elements. 
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APPENDIX A 
A. Distributed Mass Transfer Matrix for Circular Arcs 
A.1 General Remarks 
Although it is not the purpose of this thesis to derive transfer matrixes, 
the fact that the transfer matrix for a distributed mass model of a circular 
arc has not been developed and utilized before, warrants a discussion of the 
technique employed. Ref. 1 indicates the general method.* 
Since the development which follows is sae ind to a single planar 
element, it is both possible and convenient to consider an In-Plane case 
and an Out-of-Plane case. Consider the In-Plane case. A satisfactory state 
vector, consisting of three deflections and three forces, may be defined as 
follows using the notation of Sect. 2.2. 
svl 
sv2 
[SV], = | sv6 A-l 
sv12 
sv/ 
sv8 
Given such a state vector, the following equation, as suggested by 


Ref. 1, may be written: 


d(SV]ip = fA]. [sv]_, A-2 
ds 


where s indicates the length of the arc. The matrix [A] should not be 


confused with the topological matrices previously developed. 


*kcf. page 145, para. 3 of Ref. 1. 
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Matrix, [A], is an array of coefficients of six simultaneous differen- 


tial equations. It is derived in Ref. 1 and is repeated below: 


0 Dk 0 0 0 Jae 
EA 
K 0 1 0 0 0 
0 0 0 1 0 0 A-3 
[a]=|0 0 ui2 w? 02 -1 0 
0 -uw~ 0 0 0 K 
i 0 0 0 -K 0 


The parameters of the non-zero elements of [A] are defined as follows: 
K , curvature, or inverse of radius 
E, Young’s modulus 
Iz , moment of inertia about the z axis 
u, mass per unit length 
iz » radius of gyration of cross sectional area about the 
Z axis 
w, circular frequency of free vibration 
The elements of [A], for convenience, are renamed as follows after multi- 


plication by s,, total length of arc. 





Ks = 8B 

"6 a C 

EA 

Ss = H . 
oO 

EI 


vis w s* J 
~uw's =L 
[A]s. may be reordered in the following steps: 
(a) interchange rows 2 and 5 
(b) interchange columns 2 and 5 
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The following is thus obtained: 


0 0 0 0 -B C 
0 0 0 0 L B 
0 0 0 H 0 0 
[A] =0 -s J 0 0 0 A= 
S5 Oo 
B 0 Ss 0 0 0 
Oo 
Ee -B 0 0 0 0 


As indicated by Ref. 1, the eigenvalues, > , of the following array 


must now be obtained. 


[21a - [Al s)] = [sD] 


It is apparent that [SD] has the following form: 


» 0 0 +B -C 
0 a» 0 0 -L -B 
0 0 » -H 0 0 
(sD) =| ———-——- a A-5 
0 +s = ~ 0 0 
B 0 -§s 0 ~ 0 
Oo 
L +B 0 0 0 ie as 


The characteristic equation of [SD] may be obtained in a number of ways. 
The following matrix reduction scheme is most convenient. Let [SD] be parti- 


tioned into four three by three matrices as indicated above and below. 


[SD] = BE! A-6 
x | s 


of 





A reduced matrix, [SDR], having the same characteristic values, is 
given as follows: 


[sDR] = [P] - (Q]) . (s}7) . [RI A-7 


[sDR]A =~ [1] - [Q] . [R] 


The product, [Q] . [R], is given below: 


(-B2 4 cL) (-CB) — (-B) s, 
[Q] . (R] =|(2LB) (-B*) (1), A-8 
0 (-H)s (HJ) 
o 2 


Now, the elements of [SDR] } have the following form. 


On ~t +87) (cB) (B)s_ 
(-21B) (rN + B)(-L)s, A-9 
[SDR ]A = 0 (H)s| (A 7-H) 


Expansion by cofactors reveals the following polynominal as the charac- 


teristic equation, or C. E. 


2 Z 


CoE. = ( X~ + B-CL)( A74B?)( \7-HI)4( A748 -CL) (LH) 8. + 
(2LB) (CB) ( \7-HJ) - (2LB) (H) (B)s~ seins 


Further expansion reveals eq. A-l0 as the following cubic in va 


; b 
C.E. = (A‘)?-+(2B7-H5-CL) (> 7) 24(B2-2HIB-4HICLACLB “+ 


Ls”) ( \*) -(LHB“s“+HJCLB-4+CL-Hs-+HJB*) A-10a 


Further algebraic manipulation of the characteristic equation would be 
very tedious and difficult. However, numerical solution is quite simple. 
The three coefficients of the powers BES may be evaluated for an assumed 
value of frequency. The trigonometric solution of a cubic equation will 


2 
yield three values for A. The square roots of these values may then be 
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obtained yielding six values of A , generally complex. 


The solution of eq. A-2 is of the following form. 


[sv], = Saha ee [sv]. A-11 
or 
[sv], = [2], [SV]. A-lla 
where 
2), 2 Am, A-12 


Now the Cayley-Hamilton theorem may be employed as suggested in Ref. 
1. Since [A]s. is a sixth-order square matrix, the (in-plane) transfer 
matrix may be written as follows: 


* lAls 


(aig o=¢ [I] + ¢,[Als, + c,([Als,)” + c,({Als,)> + 


4 bs) 
c, (fAls 9) a cg ({A]s_) A-13 
Since the eigenvalues of [A]s_ must satisfy eq. A-13 also, six 
simultaneous equations may be written in a form similar to eq. A-13. These 


six equations may be written in matrix form as indicated below: 


a f Oo, x; BY ni ~ Co 
eX Pr WA OS Ad *, cy 
eas 1 3 x3 r3 ns X; Cy 
= < A-14 
MM) fy ry xy Ay NY AY 3 
eS Meagan Novato. We Pe c, 
ob { Ng Né ee nN x, Ce 
J 4 





or 
[ECS] = [RIA ].[C] A-1l4a 

A given value of ) _ may be represented as a complex number, R()A ) + 
JIC A ), where RC A ) indicates the real part of A and I( ZD ) indicates 
the imaginary part of » . A power of > may be indicated in a similar 
fashion as R( A - + jIC A ye When a value of A has a non-zero imaginary 
part, there will be another value of A +#which will be the complex conjugate. 
When I( A ) is zero, the corresponding rows of [ECS] and [RIA] will remain 
as shown above. 

If I( A ) is non-zero, then the two rows of [ECS] and [RIA] corres- 
ponding to the conjugate pair involved may be combined as follows to yield 
equations containing real numbers only. 


okt A D Beri Dp) fo + c,R¢ nN a =/c4R( ~~ 5 + 


e,RC A") + ccRC A”) A-15 


eC A ) sint( nN) * Co + c, ICA ) + co1( ~ 3) + c41( Dr = + 
ec, A") + eT A 7) A-15a 


These two equations may be substituted in place of the rows from which 
they were obtained so as to obtain real values only in eq. A-14 and eq. A-l4a. 
Then the coefficients may be found as 

[c] = [RIA ]7*. [ECs] A-16 

Having thus obtained the six coefficients of eq. A-13, [Z],p» may 
readily be evaluated for the assumed frequency. A similar procedure will 
yield the OQut- of-Plane transfer matrix [Z] yp: [Z] 5» and [Z])5 may then be 
combined into one twelve by twelve array which, when reordered, is the 


general three dimensional transfer matrix desired. 
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APPENDIX B 
B. Subsystem Topological Technique 
Extremely large piping systems would exceed the storage capacity of 
the large computers available. A sub-system technique is presented in this 


appendix, which may ease computer storage problems. 





Fig. B-l 


The system represented in Fig. B-l might be viewed as the four elements 
connecting nodes A and a, B and b, C and c, c and d plus two sub-systems, one 
connecting nodes a, b, c, and containing nodes 1, 2, 3, and 4, and a second 
sub-system connecting nodes d and e and containing nodes 5, 6, 7, 8, and 9. 

Let the first sub-system be considered isolated from the general system 


at nodes a, b, and c, as represented in Fig. B-2. 


-_. | 


4 


Fig. B-2 
The node stiffness matrix, [SN], may be obtained in the fashion 


previously described. This matrix would be 42 x 42. The node force matrix, 
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[FN], would have, when appropriately ordered, 18 non-zero elements in the 
low numbered positions and, since there are no forces applied to the sub- 
system at nodes other than a, b, and c; the remaining 24 elements of [FN] 
would be zero. All of the elements of the node displacement matrix [DN] 
would be generally non-zero. Let [FN] be partitioned into a non-zero 
matrix, [EN], corresponding to the forces at nodes a, b, and c, and a 
zero matrix [FN]... Let [DN], ordered so that the first 18 elements corres- 
pond to the 18 elements of [FN], > be partitioned so that the first 18 
elements make up a matrix, [DN],; and the remaining elements make up a 
second matrix, [DN], - Let [SN] be reordered and partitioned so that the 


following two equations may be written: 


[FN], = (R] . [DN], + [S$] . [DN], B-l 
[FN], = (T] . [DN], + (U] . [DN], B-2 
where 
R S 
[SN] =|} - - B-3 
T U 


Since [FN], is a zero matrix, eq. B-2 may be used to define [DN], 
as a function of [DN], > [T], and (U] as follows: 
(O] = [T]) . [DN], + (U) . [DN], B-2a 
Premultiply by the inverse of [U]. 
-(u)7) . [r). [DN], = (DN], B-4 
substitute eq. B-4 into eq. B-1l for [DN]. 


[EN], = (R] . [DN], + [Ss] . Caine sete [DN],) B-5 


ae 





or 


“a [T]) . [DN,] B-5a 


[EN], = ((R] - [Ss] . [U] 
or 
[FN], = [SN], ‘ [DN, ] B-5b 


The column matrix [FN], contains the forces applied by the sub-system 


to nodes a, b, and c. Let it be written as follows: 


[FN], = |FE B-6 








FE. DEL a 
FE, = [SN], ° DEL B-5c 
FE... DEL. 


If a similar analysis were made for the second sub-system, the follow- 


ing would result. 





= [SN], 


Similar equations can be written for those elements connecting the 


sub-systems together and to the boundaries at A, B, and C. 
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The system of Fig. B-l has been thus reduced to that of Fig. B-3. 


// a, § ¢, 13 ,e 1Y 


Fig. B-3 


The reduced node incidence matrix [A] may be written as follows: 


element Nodes 

ends A B C a b c d e 
ll, 1 0 0 0 0 0 0 0 
11 0 0 0 1 0 0 0 0 
12, 0 1 0 0 0 0 0 0 
12, 0 0 0 0 1 0 0 0 
13 0 0 0 0 0 1 0 0 
13, 0 0 0 0 0 0 1 0 
14 0 0 1 0 0 0 0 0 
14, 0 0 0 0 0 0 0 1 
1, 0 0 0 1 0 0 0 0 
1 0 0 0 0 1 0 0 0 
1. 0 0 0 0 0 1 0 0 
25 0 0 0 0 0 0 1 0 
2 0 0 0 0 0 0 0 Ji 


e = | 
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The primitive stiffness matrix of the reduced system has the following 


form: 
S44 0 0 0 0 0 
0 S49 0 0 0 0 
0 0 S 0 0 0 

[SP] = 13 

0 0 0 S34 0 0 
0 0 0 0 Ss) 0 
0 0 0 0 0 S, 

Suv Si? S43? Siu? and So are 12 x 12 arrays in this particular problem 


while Sy is an 18 x 18 array. 
The reduced system node stiffness matrix may now be formed from the 


reduced node incidence matrix and the reduced primitive stiffness matrix. 
T 


The reduced system frequency determinant, |SN| , the characteristic 
frequencies and the characteristic mode shapes may now be obtained in the 
fashion previously developed. 

In the preceding developments, trivial nodes were eliminated by using 
transfer matrix techniques. In this Appendix, we have seen that there are 
other nodes, which might be referred to as "removable nodes", which may be 


eliminated by the indicated manipulation. 
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